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Abstract 

A one-dimensional system of bosons with short-range repulsion and mid-range attraction is used 
as a laboratory to explore the evolution of many-body forces by the Similarity Renormalization 
Group (SRG). The free-space SRG is implemented for few-body systems in a symmetrized harmonic 
oscillator basis using a recursive construction analogous to no-core shell model implementations. 
This approach, which can be directly generalized to three-dimensional nuclei, is fully unitary up to 
induced ^-body forces when applied with an ^-particle basis (e.g., A-body bound-state energies 
are exactly preserved). The oscillator matrix elements for a given A can then be used in larger sys- 
tems. Errors from omitted induced many-body forces show a hierarchy of decreasing contribution 
to binding energies. An analysis of individual contributions to the growth of many-body forces 
demonstrates such a hierarchy and provides an understanding of its origins. 
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I. INTRODUCTION 



A major goal of nuclear structure theory is to make quantitative calculations of nuclear 
observables starting from microscopic internucleon forces. The Similarity Renormalization 
Group (SRG) QHi, through a continuous series of unitary transformations that soften 
initial interactions, can dramatically reduce the computational requirements of low-energy 
many-body calculations [5, IB, @1- At the same time, the SRG induces many-body forces 
as it evolves the Hamiltonian. For the SRG to be a useful tool, we must develop methods 
for calculating these many-body interactions and establish the conditions under which an 
initial hierarchy of many-body forces is maintained. In this paper, we study one-dimensional 
systems of bosons as a proof-of-principle of a practical method to evolve and evaluate such 
forces, and establish a road map for full three-dimensional calculations. 

The SRG can be implemented as a flow equation for the evolving Hamiltonian Hg, 

^ = [r^s.H:\ = [[G,,H,iH,]. (1) 

Here s is a flow parameter and the flow operator Gs specifies the type of SRG 0, 0]. Most 
previous apphcations to nuclear structure have been in a momentum basis, where decoupling 
between low-energy and high-energy matrix elements is naturally achieved by choosing a 
momentum-diagonal flow operator such as the kinetic energy Trei. However we can evaluate 
Eq. ([1]) with Gg = Tj-ei in any convenient basis. 

In Ref. jsl, a diagrammatic approach to the SRG equation was introduced, which orga- 
nized the independent evolution of two- and three-body (and higher-body) potentials. This 
formalism is necessary in a momentum basis to avoid "dangerous" delta functions from spec- 
tator particles. An alternative is to work in a discrete basis for which Eq. ([T]) can be applied 
directly in each n-body sector without a separate evolution for few-body forces; that is, the 
Hamiltonian is evolved as a whole. We adopt this approach in the present work, using har- 
monic oscillator wave functions as our basis and mimicking the formalism used in the no-core 
shell model (NCSM) [q, T^, 11] to create properly symmetrized (for bosons) matrix elements 



in relative (Jacobi) coordinates. The restriction to one dimension makes the construction 
particularly straightforward and requires only moderate matrix sizes. We use bosons in this 
paper for easy comparison with existing model analysis. The boson ground states coin- 
cide with fermion ground states when the flavor degeneracy is greater than the number of 
particles, because the overall anti-symmetrization is realized by the flavor wavefunction. 

We use simple flavor-independent potentials that imitate the short-range repulsion and 
mid-range attraction characteristic of realistic local nuclear potentials. Previous studies 
of the SRG imply that properties of the transformations are primarily due to the matrix 
structure (Gs, Hg, choice of basis, etc), so we expect to be able to directly carry over 
at least some of our observations to three dimensions. Because the NCSM formalism is 
already developed, the generalization to three-dimensional fermionic calculations with spin- 
isospin degrees of freedom and using realistic nuclear interactions should be algebraically 
straightforward (although far more computationally intensive). 

The plan of the paper is as follows. In Section [Tll we develop the one-dimensional version 
of the oscillator basis techniques. These include building symmetrizer operators to construct 
symmetric basis states of bosons and fermions, embedding potential and kinetic energy 
operators in a given oscillator basis, and developing a scheme to organize the oscillator basis 
states in a universal and scalable manner. In Section llllt results are given for two- through 
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five-particle systems. The two-particle calculations establish that the models simulate the 



relevant features of previous nuclear calculations 0, H, 0, 12]. Then three- and higher- 
particle calculations are used to explore the running of the few-body potentials. Section [IV] 
summarizes our observations and conclusions. 



II. FORMALISM 



In this section, we adapt to one dimension the recursive symmetrization formalism de- 
veloped by Navratil and collaborators for use in NCSM calculations with a translationally 
invariant harmonic oscillator basis [sl, [13, 11 1. While the three-dimensional formalism is well 



documented, the one- dimensional analog is not, so we provide a self-contained treatment 
here. 



A. Jacobi Coordinates 

The initial (i.e., unevolved) one-dimensional Hamiltonian for A bosons of equal mass m 
with a local two-body potential has the first-quantized form (in units with h= 1) 

H=^Y.kl+ J2 V{x.-x,), (2) 

i=l i<j=l 

where the Xi are single-particle coordinates and the ki are single-particle momenta. To 
connect to the nuclear problem of interest that uses potentials given in a momentum basis 
(e.g., chiral effective field theory potentials), we calculate matrix elements using harmonic 
oscillator basis states in Jacobi momentum coordinates. This representation also provides a 
clean visual interpretation of the SRG evolution of potentials. 

With equal-mass particles, a convenient set of relative momentum Jacobi coordinates is 
defined by (for j = 1 to A — 1) 



Pi 




T ( 7 5Z " ^i+i I , (3) 



3 ^=l 



where the ki are the single-particle momenta of the A particles. Sometimes for convenience 
we will use p = pi and q = p2 when restricted to A < 3. We define the Fourier transform 
V{xi — X2) to V{pi,p[) as 

Vipi,p[) = J ViV2i,)e-'^P^~^'^^'^di, , (4) 

where ii = (xi — X2)/V2 is the coordinate conjugate to the Jacobi momentum pi. We 
introduce a set of harmonic oscillator states \nj) corresponding to each of the coordinates 
of Eq. ([3]), so a general product basis state has the form 

A-l 

with rij = 0, 1, 2, ... , iVmax for each j. In the next section we discuss how to build linear 
combinations of these states that have the appropriate symmetry. 
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B. Symmetrization 



We carry out the SRG evolution for each A-particle subsystem in a complete basis of 
properly symmetrized states, which will be linear combinations of the basis states of Eq. (Q. 
The symmetrization procedure is adapted from the procedure developed for NCSM calcu- 
lations 0, 11]. This entails symmetrizing first the two-particle system and then using a 



recursive procedure to go from the (A — l)-particle basis to an A-particle basis. At each stage 
we keep only symmetric states, identified as eigenstates of the symmetrizer with eigenvalue 
unity. 

The two-particle system is specified by the oscillator number ni. The symmetrizer is 
(l-|-Pi2)/2, where Pij is the exchange operator between particles i and j. Because Pi2|'^i) = 
(— 1)"^ |ni), the symmetrizer in the two-particle case has eigenvalue one acting on states with 
rii even and zero when acting on states with ni odd. Thus the symmetric basis states have 
Til even and we simply omit the odd states. Following conventions from Ref. [9|, we label 
these eigenstates as \N2i2), where N2 is the total oscillator number of the symmetric state 
and ^2 is an arbitrary label which distinguishes states degenerate in N2. In the two-particle 
case the notation is trivial, with A^2 = ni even and ^2 = 1. We write eigenstate projection 
coefficients as {N2i2\\ni) = 6N2,ni{^ + ("~l)"0/2- These are referred to in the literature as 
the "coefficients of fractional parentage" . 

A three-particle basis is specified by product states of the two-body symmetric eigen- 
states, \N2i2), and single-particle states with the oscillator number corresponding to the 
third particle, 1^2). The symmetrizer for this system is governed by the permutation group, 
5*3, which can be defined by just two of its generators. Here we choose the permutation 
operators P12 and P23. The symmetrization operator can be written as 

5 = i(l + P12 + P23 + P12P23 + P23P12 + P12P23P12) . (6) 

We build this symmetrizer in the basis |A^2^2;^2) = |A'^2^2)|^2) where the states |A^2^2) are 
already eigenstates of Pu with eigenvalue one, so Eq. ([6]) reduces to S* = (1 + 2P23)/3. 
In this basis, the matrix elements of P23 can be expressed as 

{N^i2; n2\P23\N2i2; 112) = 5N',N{n[n2\nin2)3 , (7) 

where N = N2 + n2 = A^2 + ^2 ~ ^' (^i'^2 1^1^2)3 is the one-dimensional harmonic 
oscillator transformation bracket for particles with mass ratio 3. We construct these trans- 
formation brackets and generalize to mass ratio d in Appendix |Al By diagonalizing this 
symmetrizer we identify the symmetric eigenstates of the system as the ones with eigenvalue 
unity. We keep only those states and discard the others. This set of eigenvectors gives us 
the coefficients of fractional parentage, {N2i2] n2\\N3i3), of the three-boson symmetric eigen- 
states, \N3i3), in terms of the original partially symmetrized three-particle space, |A^2^2; ^2). 
Note that is not trivial like 12, because in the three-body system there are eigenstates 
degenerate in the total oscillator number, N2 + n2. The label keeps track of those de- 
generacies. We find in the one-dimensional system of bosons that the fraction of symmetric 
basis states for A = 3 is about one-fifth. For A = 4 the reduction in number of states is 
above 90%. 

To construct the basis states for higher A, we generalize this procedure. To go from 
A — 1 to A we need only to symmetrize between the last two particles, so we construct the 
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symmetrizer 

'^a = ^(1 + (A-1)P(a-i)a) (8) 

in the space of {A — l)-particle symmetric eigenstates and the additional Jacobi state, ua-i- 
We label the basis of this space as |A^a-iM-i; ^A-i)- The matrix element of the exchange 
operator in this space is 

{N'^-ii'A-v n'A-i\P{A-i)A\NA-i'iA-i; nA~i) 

X (A^A-2M-2;?^A-2||A^A-1M-i)(^A-2^A-iI'^A-2^A-i)a(A-2) , (9) 

where the sum is over Na-2, 'iA-2, nA-2 and n\_2. The only significant difference from 
the three-particle case is that we must sum over the components of the A — 1 subcluster 
symmetric states to get all the contributions to the exchange of the last two bosons, nA-2 and 
ha-i- The parameter d = A{A — 2) can be derived by taking the last two Jacobi coordinates 
Pa-2 and Pa-i, as defined in Eq. ([3]), and finding the transformation that exchanges particles 
labeled by kA-i and kA- This procedure is shown in Appendix lAl 

For fermions, we need a complete basis of fully anti-symmetrized states. If we consider 
the one-flavor case, the procedure for our one- dimensional model is a trivial modification of 
Eqs.(in]) and ([S]), namely all odd permutations come with a minus sign. Thus, for A = 3 the 
anti-symmetrizer can be written 

A=^(1-2P23), (10) 

where P23 acts on the flavor space as well. If there are more flavors than particles and 
the interaction is flavor independent, the spatial wavefunction for the ground state will be 
symmetric and correspond to our boson ground state wavefunctions. For realistic three- 
dimensional nuclei, the required construction of an anti-symmetric Jacobi basis with full 
angular momentum coupling has been worked out for the NCSM by Navratil et al. jo], [l^, 11 



C. Hamiltonian Matrix Elements 

To obtain the Hamiltonian in the symmetric eigenbasis of the general A-particle system, 
we employ a recursive embedding procedure that utilizes the partially symmetric bases 
developed for the symmetrization operator. First we treat the kinetic energy and then the 
potential. 

The relative kinetic energy in the three-particle system is the total minus the center-of- 
mass kinetic energies: 

(fci + k2 + ksf 
2(3m) 

(11) 

where the pj's are defined in Eq. ([3]). Momentum basis states are organized by increasing 
kinetic energy. We can project Tj-ei directly onto the three-particle oscillator basis by using 
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the ladder operator definitions of tlie Jacobi momenta. Tlie projection of T^^i into the \n1n2) 
basis is 

I ! ! \rr \ \ / / / iPl + P2 I \ 

{n^n2\Tre\\nin2) = {n^n^l \n1n2} 

= ^-y^^^^aK^I - ^1)' + (^2 - V2y\nin2) , (12) 
where the rji and ri2 operators act on the rii and 77,2 spaces, respectively. Continuing, we get 

(n>2|Trei|nin2) = 2L ^^ IKK^i - Vi)^\ni)Sn^,n'^ + {n2\{vl - V2y\n2)Sn,,n'^] 

- {2ni + l)6n[,ni) ^n'2,n2 

+ (a/ (^2 + l)(n2 + 2) 5„'2,„2+2 + ^n2{n2 - l)5^,„2-2 

- (2^2 + l)5n^,n2) ^ni,ni • (13) 

As noted, we keep only the ni-even states using the projector {N2i2\n\), and we can sym- 
metrize the full three-particle system with the symmetric eigenstates, jiVsis) whose compo- 
nents are given by {Nj,%j,\N2%2\ ^2). 

To derive the A-body kinetic energy in the symmetrized basis, (Tyi)syni, we use a recursive 
procedure on the {A — l)-body result to find the A-particle space operator matrix elements: 

A-l 

1=1 

= {K\\K-inU-i){K-inU~i\TA\NA-inA~i){NA^inA-i\\NA) 
= {N'A\\N'A_,n'A^,)[{N'A_,n'A.,\{TA-i)^^^ 

= (^^II^A~l^A-l) [{NA^l\{TA-l)sym\NA-l)Sn[^_^,nA-i 

+ ^Jv;i_i,iVA-i("'A-ibi-i/2"^l'^A-i)] {NA~inA-i\\NA) 



-6„ 



-2,nA- 



^1 - Sn'^_^-2,nA-i Vi^A-i + l)(n^-i + 2) )] 



X {Na-iUa-iWNa) , (14) 

where we have suppressed the m's and i^-i's for simplicity after the first line. Intermediate 
summations over Na-i, tia-i, «a-i, N'a-h ''^'a-i; ^a-i implicit. 

In the same manner as the kinetic energy, we can recursively embed the potential in 
the A-particle space, starting with the two-body interaction between the first two particles. 
Because we are working in fully symmetrized few-body spaces we do not need to consider all 
pair-wise interactions, but only one such pair and scale by the number of interactions. For 
instance, in the three-particle system the full two-body interaction is V^^^ = V12 + V23 + V13 = 
3V^i2. In a general A-particle space, this becomes V^^^-* = (^)Vi2. The matrix element 
of a two-body potential, V12, in the relative coordinate harmonic oscillator basis, |?t,i), is 
(n'j^l Vi2|ni) = J {n'^\p'i) {p' \V12\p) {p\ni) dp dp' where the matrix elements of {p'\Vi2\p) are given 
by Eq. dD. 
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Once in the oscillator basis, embedding in a larger particle space is a straightforward pro- 
cess. Starting with the two-body interaction, V12, the two-body oscillator symmetric states 
are isolated using the projector (iV222|^i) which picks out just the ni-even states. Embed- 
ding this interaction in the three-particle space involves adding a new Jacobi coordinate, 
|n2), to the existing system. With respect to the two-body interaction, V12, this additional 
coordinate is associated with a delta function, 6n2,n'^- Finally we obtain the symmetric three- 
particle states by using the projector, (A^aisl |Ai'2«2; ^2)- Multiplying by (2) = 3 gives us the 
full strength of the two-body interaction. 

In general we can write this procedure as an expansion of the final A-particle symmetric 
space matrix elements of V12: 

= (iV;i||iV;_in:4_i)(iV;,_y^_i|VA|iV^-inA-i)(iV^-inA-i||iV^) 

X {N'A_^nj^_^\{VA-l)sju,5n'^^^,nA^i\^A-inA~l){NA-inA~l\\NA) (15) 

where again we have dropped the Zyi's after the first line for simplicity and intermediate sums 
are implicit. We remind the reader that ua-i can only take values from to N — Na-i-, where 
is the total oscillator quantum number used to organize the states. We start with the 
two-particle space and work our way up to the A-body space, embedding the interactions 
successively in each sector using Eq. fllSp . When symmetrizing Va we must embed the 
symmetrized Va-i with the appropriate combinatoric factor included as explained above. 
This factor derives from the fact that we had embedded a 2-body force in the A — 1 space 
that is now to be extracted and embedded in the A-particle space. Thus we must remove the 
old factor ("^2 ^) and multiply by the new (^) factor, which has the net effect of multiplying 
hYA/{A-2). 

Any initial three-body force (discussed below) is embedded in the same manner as above 
except that it originates in the three-particle space. The initial three-body force is a function 
of two Jacobi momenta, which we transform directly into the partially symmetrized three- 
particle oscillator space and then use all of the same embedding procedures developed above. 
Note that two- and three-body forces must be embedded in higher spaces with different 
symmetry factors, (^) and (3) respectively. 

In previous formulations of this recursive approach [9|, subsequent potential embeddings 
are achieved by making a change of coordinates for the last two Jacobi momenta. For 
systems with A > 5, the three-body force requires a similar change of coordinates for the 
last three Jacobi momenta. Such a scheme is unnecessary here. 

D. SRG Evolution 

Once we have constructed a complete symmetrized basis for A particles (specified by the 
value of A^max) and evaluated the Hamiltonian matrix elements in this basis, applying Eq. ([1]) 
with Gs = Trei is immediate. That is, we have coupled, first-order differential equations for 
each matrix element of the Hamiltonian, with the right side of each equation given by a 
series of matrix multiplications. This is efficiently implemented in any computer language 
with matrix operations and differential equation solvers. 
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Individual matrix elements of the Hamiltonian obey the SRG's differential equations: 

^JN'^'^\{VA)s\NA^A) = {N'Ai'A\[[T,H,],H,]\NAtA) 

= {N'xi'a\TH,H,\NaIa) + {N'xi'a\HsHsT\NaIa) 

-2{N'xi'a\HsTHs\NaIa) . (16) 

We have defined dT/ds = so that all of the flow occurs in the matrix representation of the 
potential, {Va)s- Using the matrix representations of T and Hg in the {NaIa) basis, the right 
side of Eg. (1161) is simply a series of matrix multiplications. The initial condition at s = 
is the initial Hamiltonian, {N'^i'^\T + Va|A^a^^), which can have few-body components in 
Va- We consider here both two-body-only and two-body plus a three-body component (for 
A>3). In practice we often use the flow variable A = because there is no explicit s 

dependence in Eq. ([T]), switching variables is trivial. 

To carry out the SRG evolution, we use a built-in MATLAB differential equation solver, 
such as the MATLAB function ode23, which is an implementation of a Runge-Kutta differ- 
ential equation algorithm. We studied the running time to evolve the various SRG schemes 
(i.e., the choice of Gg) by plotting the time to run versus the evolution parameter, s. We 
find a straight line as s increases, indicating no stiffness, in every combination of potential 
and SRG scheme used to date. 

The SRG induces few-body forces as it evolves an initial interaction in a few-particle space. 
To study the contributions of different few-body forces we must isolate these components 
of the full interaction. By definition the two-body force evolution keeps the A = 2 binding 
energy invariant under evolution. We can isolate the three-body force from the two-body 
matrix elements by embedding the evolved two-body-only force in the three-particle space 
and subtracting it from the full two-plus-three-body evolved interaction. In our MATLAB 
implementation these procedures take only a few lines of code. 



III. RESULTS 

A. Initial ("Bare") Interactions 

The bulk of our calculations adopt a model from Ref. [is!] that uses a sum of two gaus- 
sians to simulate repulsive short-range and attractive mid-range nucleon-nucleon two-body 
potentials: 

V(^){x) = ^e"-^/'^? + ^e-^^/'^i (17) 

or 

V^'\p,p') = J^e-(P-^')''^?/« + J^e-(^-P')'-^'/« . (18) 
27rv2 27rv2 



The parameters used in Ref. were chosen so that the one-dimensional saturation prop- 
erties correspond to empirical three-dimensional properties, but we also want to explore a 
range of parameters to test what behavior is general and what relies on specific features. 
We start with the parameters listed in Table [H The potential Va is from Ref. and is 
plotted in Fig. [TJ We will fix the range of the attractive part and vary the relative strength 
and range of the repulsive parts and visa versa. We also vary the purely attractive poten- 
tial Vp, which was used in Ref. [ij] and is also plotted in Fig. [H The eigenvalue problem 
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TABLE I: Parameter sets for the two-body potential of Eq. (jlSp . 



name 


Vi 


0-1 


V2 


0-2 




12. 


0.2 


-12. 


0.8 




0. 


0.0 


-2.0 


0.8 



for the relatively small matrices considered here can be solved by any conventional matrix 
diagonalization program (MATLAB was used here). 
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FIG. 1: (color online) Potentials (dotted line, with axis on right) and probability distributions 
(other lines, with axis on left) for the lowest two-body bound state as a function of x = — a;2| 
at different stages in the SRG evolution (A = The left plot is Va and the right plot is Vg 

(see Table H). 



To test that the symmetrized harmonic oscillator basis was correctly constructed for 
A = 2, 3, and 4 (see below), we first diagonalized the Hamiltonian using the purely attractive 
gaussian two-body potential Vg. The normalization is such that Vp{x) becomes a delta 
function with strength V2 as ^ [3] (note the numerical factors from the Fourier 
transform because of our normalization of the Jacobi momenta). This limiting case has 
a known analytic solution for the (only) bound state of A bosons. For finite (T2, we were 
able to confirm the accuracy of the diagonalizations as a function of the basis size A^max by 
comparison to coordinate-space stochastic variational method (SVM) calculations using a 
published code adapted to one dimension 16 . 

We also explore the evolution of Hamiltonians with an initial three-body force. We choose 
a regulated contact interaction in the three-particle momentum space. 



V^^\p, q,p\ q) = CEfAip, q)fk{p\ q) 



(19) 



where ce is the strength of the interaction and the form factor /a depends on the Jacobi 
momenta as 

-{(p2+92)/A2)" 



/A(p,g) 



(20) 
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TABLE II: Ground-state energies for two-body potentials from Table |I] with various strengths of 
the initial three-body potential Eqs. ([T9|)-(l20]) with A = 2 and n = 4 for ^4 = 2, 3, and 4. 



y(2) 


CE 


E2 


E3 


Ea 


Va 


-0.10 


-0.920 


-3.223 


-7.125 




-0.05 


-0.920 


-2.884 


-5.832 


Vo. 


-0.01 


-0.920 


-2.628 


-4.906 


Va 


0.00 


-0.920 


-2.567 


-4.695 


Va 


0.01 


-0.920 


-2.507 


-4.494 




0.05 


-0.920 


-2.278 


-3.798 




0.10 


-0.920 


-2.027 


-3.179 


Vp 


-0.10 


-0.474 


-3.379 


-8.412 




-0.05 


-0.474 


-2.283 


-5.727 


yp 


-0.01 


-0.474 


-1.792 


-4.183 


Vp 


0.00 


-0.474 


-1.708 


-3.846 




0.01 


-0.474 


-1.626 


-3.517 


Vp 


0.05 


-0.474 


-1.370 


-2.451 




0.10 


-0.474 


-1.240 


-1.874 



The regulator cutoff A sets the scale of the fall-off in momentum and n determines the 
sharpness of this fall-off. This form is analogous to the regulated three-body contact inter- 
actions used in chiral effective field theory [17]. We have not explored in detail the impact 
of adjusting A and n but have focused on how the SRG handles a varying strength ce- All 
results here are for A = 2 and n = 4. 

Most of the figures in this paper show calculations with A^max = 28. With this basis size, 
the ground-state energies are generally converged to one part in 10"^, which is more than 
sufficient for our purposes. As usual, increasing A^max leads to rapidly increasing matrix sizes 
and computation times; times for A = 3 with A^max = 32 are a factor of 3 longer than with 
A'max = 28 and with A^max = 40 the time increases by another factor of 10. A sampling of 
ground-state energies are given in Table [Tll 

B. Two-body Results 

We first consider the bound state of two identical bosons using the potential V^. Because 
the SRG is a series of unitary transformations, we expect that the binding energy will not 
be changed by evolving the two-body interaction in the two-particle space. Indeed, we 
find it to be constant to high accuracy. The ground-state wave function, however, changes 
dramatically, as seen from the probability densities plotted in Fig. [H The initial probability 
density exhibits a sizable "wound" near the origin that is filled in as A decreases. By A = 2 
there is no signature of a repulsive core (and the wave function is modified out to larger x) . 
This is the same pattern seen for the S-wave component of deuteron wave functions starting 
from three-dimensional nucleon-nucleon S-wave potentials with strong repulsive cores such 
as Argonne fig, with the "uncorrelated" final wave function at A = 2 roughly comparable 
to A = 1.5 fm~^ for the deuteron [5|. 

The evolution of the potential in the momentum basis, shown as a color contour plot in 
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FIG. 2: (color online) Even part of the SRG potential [V^ {p,p') + Vs {p, —p')] in dimensionless 
units as a function of p and p' for A = oo, 5, 3, and 2 (where A = The initial potential is 

Va from Table HI 

Fig. [2], also demonstrates this behavior. (The even part of the potential is shown, which 
is the analog of the S-wave part.) The initial potential is dominated by strongly repulsive 
matrix elements coupling low and high momenta. The evolution in A band diagonalizes the 
potential to a width in of roughly A^ while a soft attractive part emerges in the low- 
momentum region. The pattern in Fig. [2] reflects increasing non-locality as A is lowered, 
which in turn reduces the wound in the wave function. From the probability density and 
the momentum space plots we estimate that evolving to halfway between A = 2 and 3 for Va 
corresponds roughly to the A scale typically used in nuclear structure calculations (around 
2fm-^). 

C. Three-body Results 

To calculate properties of the three-particle system we construct the Hamiltonian in the 
basis of symmetric three-particle eigenstates as described in Sec. Ill CI The SRG evolution of 
the potential in the three-particle space leaves the ground state energy invariant if the full 
Hamiltonian is kept, because the transformations are unitary. However, the Hamiltonian 
matrix elements in this space do not follow simply from the pairwise sum of the two-body 
potential matrix elements; as the SRG evolves a three-body force is induced even if its initial 
strength is zero. 

The effect of this full three-particle space SRG evolution is shown in Fig. [3] for initial 
two-body potential Va and with initial V^*^^-* = (non-zero values of ce are considered in 
the next section). We plot the ground-state energy for the three-particle system both with 
the initial two-body interaction embedded in the three-particle symmetric space and then 
evolved (the red curve with squares) and also with the two-body interaction evolved in the 
two-particle space before embedding in the three-particle space at each A (the black curve 
with circles). We can see that the energy evaluated with the two-body interaction alone 
deviates noticeably as A drops below 10. This variation is the signature that the two-body 
transformation is only approximately unitary in the three-particle sector. The error reaches 
a peak in A between 2 and 3 and then decreases. The same pattern has been observed for 
NN potentials in three dimensions |6|] and remains qualitatively the same when parameters 
in the potential are varied (e.g., see the right plot in Fig. [3]). We made the same calculation 
using the purely attractive initial two-body potential Vg, which is shown in Fig. HI Here 
the induced three-body force has the opposite sign and there is no maximum, which implies 
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FIG. 3: (color online) The lowest bound-state energy £"3 for a three-particle system as a function 
of A with the initial two-body-only potential Va- The (red) curves with squares include the full 
evolution of the Hamiltonian while the (black) curves with circles use the two-body potential 
evolved in the two-particle system. The right frame shows two additional results from varying cji 
and V2 from the values in Table HI 



Orr 




FIG. 4: (color online) Same as Fig. [3] but with initial potential Vg. The right frame shows two 
additional results from varying cj2 and V2 from the values in Table HI 

that the qualitative pattern of evolution is dictated by the interplay between attractive long- 
range and repulsive short-range parts of the potential. These features are explored further 
in Sect. HTTEI 

In Fig. [5] we test SRG decoupling [12] within the harmonic oscillator basis. An initial 
two-body potential is evolved in a large A = 3 space both with and without the induced 
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10 20 30 40 '0 10 20 30 40 



N , N , 

cut cut 

FIG. 5: (color online) Decoupling in the three-particle system. The intial Va potential is evolved to 
each A shown in a basis with A'^max = 40. On the left only the two-body potential is kept while the 
full potential is kept on the right. Matrix elements of the potential are set to zero if one or both 
states have > A'cut and the resulting Hamiltonian is diagonalized to obtain the ground-state 
energies plotted. 

three-body interaction. Then the Hamiltonians for selected A values are diagonalized in 
bases of decreasing size, as measured by 'Wcut" , which is the cut-off applied to the potential 
to study its decoupling properties. (That is, the potential is set to zero for matrix elements 
for which one or both states has > Nc^f) The left panel shows the results when only the 
two-body evolved potential is used. The degree of decoupling is measured by the point of 
departure from the asymptotic energy for A/'max = 40 as A"cut is lowered. As the potential 
is evolved from the initial potential (A = oo) down to A = 2.5, decoupling is achieved for 
smaller spaces, which means convergence is reached for smaller basis sizes. This is the same 
pattern as found for realistic NN potentials [l^ . 

In the left panel of Fig. [5l the ground-state energies asymptote to different values because 
the evolution is not completely unitary. One might imagine that this would affect the 
decoupling, but we show this is not the case here. In the right panel, the induced three- 
body interaction is kept, so the curves asymptote to the same energy at large A^max, while 
the same pattern of decoupling is observed. We note that the decoupling benefits afforded 
by evolution in the oscillator basis are less straightforward than in the two-body momentum 
basis studied for NN in In particular the cut-off errors increase for A smaller than 

the point at which the two-body-only binding energy is at a minimum (i.e., for A = 1.5 in 

Fig. ED. 

D. Results for ^ = 4 and A = 5 

Next we turn to A = A and A = 5, where we expect to see the effects of induced three-, 
four- and five-body forces. The key issue is the relative sizes of these contributions; we 
are looking to test whether an initial hierarchy of few-body interactions is preserved and 
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FIG. 6: (color online) The lowest bound-state energy E4 for a four-particle system as a function of A 
with the inital two-body potential Va and different initial three-body force strengths (cg = ±0.05). 



therefore can be truncated with a controlled error. (Of course, generalizing these results to 
three dimensions will require repeating these tests for real 3D nuclei.) 

In applying the SRG in the four-particle space we have three different calculations of the 
ground-state energy to compare. The first is the two-body potential embedded successively 
in the three- and four-particle spaces and then evolved in the four-particle space. The 
resulting unitary transformations will induce three- and four-body interactions that leave 
the eigenvalues invariant. We can also evolve in the two-particle space before embedding 
in the three- and four-particle spaces and diagonalizing. As we saw before in Fig. [3] and 
see now in Fig. [6l the two-body-only evolution deviates because the Hamiltonian is not 
evolved by an exactly unitary transformation. Finally, to find the relative size of the three 
and four-body interactions we can evolve in the three-particle space, thereby inducing only 
three-body forces. Note that the two and three-body forces must be embedded differently 
in the four particle space because they have different combinatoric factors associated with 
them, i.e., there are (2) = 6 pairs and (3) = 4 triplets. So, the proper mixture of two- and 
three-body force contributions to the four-particle system interaction is = GV^"^^ + 4^^*^^). 

All three of these calculations for A = 4 are shown in Fig. [HI for the two-body potential Va 
and several choices of the initial three-body force. The magnitude of ce was chosen so that 
the fractions of the A = 3 and A = 4 ground-state energies from the three-body interaction 
are roughly comparable to the corresponding fractions for nuclei using typical realistic NN 
potentials. The qualitative behavior is similar for other choices of ce and V'^'^\ In all plots 
the curves for the two-body-only (black line with circles), the two-plus-three (red line with 
squares), and the full two-plus-three-plus- four interaction (blue line with diamonds) show the 
hierarchy of different few-body interaction components. The difference between the square 
and diamond lines represents the contribution of the four-body force, and the difference 
between the circle and square lines is the contribution of the three-body force alone. 

The four-body contribution is at most ten percent that of the three-body, which is itself 
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FIG. 7: (color online) The lowest bound-state energy for a 5-particle system as a function of A 
with an initial two-body-only Va potential for several values of A^max- 

small compared to the two-body contribution except when the latter gets small for small A 
(note the expanded scales on the figures). Considering calculations with different ce values, 
we see that the A dependence of the induced four-body part depends on the interplay of 
initial and induced forces. In some cases noticeable (but small) evolution starts at A = 10 
while in other cases it is deferred until much smaller A. Regardless of the details, we stress 
that there is no sign that induced many-body forces have rapid growth with A or exhibit 
unusual scaling. 

We repeated for A = 4 our test of decoupling that was shown in Fig. O for A = 3. A 
similar pattern of decoupling is found, namely an increased degree of decoupling until a 
A corresponding to the minimum of the two-body-only ground-state energy of the A = 4 
system, after which it deteriorates. 

In Fig. [7] we show results for the SRG evolution, with initial potential Va and no initial 
few-body interactions, in a five-particle system for several values of A^max- The right panel, 
with A'max = 28, shows the converged result. We see a decreasing hierarchy of induced 
many-body contributions for all A; the five-body contribution is essentially negligible (or 
not distinguishable from numerical noise). Differences at the lower Amax sized spaces arise 
both because the space needs to be large enough for convergence to the exact energy eigen- 
values but also because the initial evolution of the potential needs a sufficiently large space. 
Decoupling may improve this feature but is dependent on the type of SRG used and the 
basis in which it is implemented. This is important to keep in mind for the generalization 
to the realistic nuclear problem. 



E. Analysis of Three-Body Force Running 

In this section we examine how induced three-body interactions evolve in our one- 
dimensional laboratory and make connections to the diagrammatic expansion developed 
in Ref . jl] . As raw material for this analysis, we plot in Fig. [8] the error of the evolving 
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FIG. 8: (color online) Differences of two-body-only and two-plus-three-body A = 3 ground-state 
energies as a function of A. Each of the parameters of the potential Va are varied in each plot as 
the other parameters are held constant. The upper panels vary the ranges while the lower vary the 
strengths; the left vary the attractive part and the right vary the repulsive part. 
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FIG. 9: (color online) Same as Fig. [8] but for Vg 
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FIG. 10: A diagrammatic decomposition of the SRG equation (j2ip . A circle at a vertex denotes a 
commutator with Tj-ei. 

two-body-only binding energy while varying several parameters of the initial two-body in- 
teraction Va for A = 3 ground-state energies. In the left plots we vary the range (top) and 
the strength (bottom) of the attractive part of V^. In the right plots we vary the range 
(top) and the strength (bottom) of the repulsive part. In Fig. [9] we present similar plots for 
a simpler system that starts with the initial attraction-only potential Vg. Note that what is 
plotted shows the evolution with A of the ground-state expectation value of the three-body 
force. 

Certain qualitative features are found as expected in these figures. Shorter ranges imply 
enhanced coupling from low to high momentum and therefore we anticipate that the evolu- 
tion will start sooner (i.e., at higher A). This is seen clearly on the left in Fig. [9] and for the 
variation of the (shorter-ranged) repulsive potential in Fig. [8] (top right plot). There is also 
an unsurprising increase in the magnitude of the induced three-body interaction at each A 
with increased magnitude of the potential, as seen on the right in Fig. [9] and on the bottom 
left in Fig. [HI For a more definitive analysis we need to recall the discussion from Ref . . 

The SRG evolution equation for the three-particle sector in the notation of Ref. is 

where each bar denotes a commutator with Trd. We remind the reader that dT^-d/ds = 
by construction. A diagrammatic decomposition of this equation is shown in Fig. [TDl In 
the two-body sector, the equation reduces to the first term on the left and the first and 
third terms on the right (the first row in Fig. [TU]) . These terms keep two-particle energy 
eigenvalues invariant under evolution. In the three-particle sector, Eq. (1211) results in not 
only these two-body graphs with a disconnected spectator but additional graphs involving 
connected combinations of two and three-body interactions. The diagrams with two-body 
interactions and a disconnected spectator line satisfy the two-body evolution equations, and 
so will cancel out of the full three-particle-sector evolution equation. Thus the evolution of 
the three-body interaction is dictated by the connected diagrams (the second row in Fig. [TOl) . 
In summary, the evolution of the A-hodj potential in the A-particle system is given by 



dV^ 
ds 



[Vs,H,]a, (22) 



where the "A" subscript on the right side means the fully connected A-particle terms. 

To make a connection between the individual terms in the three-body interaction evolu- 
tion and the running of the ground-state energy, we need to derive the evolution equations for 
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(3) 

the expectation value of Vs in the ground state. Denoting the ground-state wave function 
for the A-particle system by \ip^), it evolves according to 

= Us\^|Jt,) , ^J^Pf)=Vs\^Pt), (23) 
where Us is the SRG unitary transformation at s and 

Vs = = -vl ■ (24) 

Then the matrix element of an operator Og evolves according to 

i(^s\Os\i^s) = ii^fl^ - [Vs, . (25) 

If the operator Og evolves according to Os = UsOs=oUl, then the matrix element vanishes, 
as when Og = Hg. 

However, if we wish to see how one part of Hg evolves, such as the expectation value of 
V^^\ we obtain 

= - [Vs, Vj^i^t) , (26) 

which does not give zero in general because Vj^^ ^ UgVgtlWg. In the two-particle case, the 
analog of Eq. fl26l) gives d{V^'^^)/ds = ([r/s, Trd]). In the three-particle case, we can expand 
Eq. (EE]) as 



^Ji^sivj'^li^g) = mVs,Hg]s-[Vs,V}'^]\ijg) 

= mvf\Ze^] + [vf\vj% + [vf\vj'^] + [vf\vj^^] + [vf\vj'^] 

-[vf\vj'^]-[vf\v}'^m 
= mv?\Hg] + [vT,v}\ - [vf\vjni^g) 

= mvf\Vg^\-[vf\vj''>Mg), (27) 

where vf^ and V^^^'* are the commutators vf^ = [Trd, V^i^''] and vf^ = [T,.ei,vj^'']. In the 

third line, the expectation value of the commutator, [V^^ ,Hg], vanishes identically. 

The term with the subscript "c" has had the two-body disconnected diagrams removed. 
In our MATLAB implementation, this subtraction is achieved by first embedding the two- 
particle-space evolved version of this commutator in the three-particle space. Computing 

[Vg , V^'-^''] in the two-particle space alone involves only the one-loop two-body interactions, 
so embedding in the three-particle sector results in only the disconnected parts. This dis- 
connected part can then be subtracted from the total three-particle sector version of the 
same commutator, leaving only the three-particle fully connected part. 

It is most useful for our analysis to convert from derivatives with respect to s to derivatives 
with respect to A using ^ = -^^S- HH we show the ground-state expectation 

values of the right side of Eq. fl26|) . which are broken down into the two terms from the 
right side of Eq. fl27j) for A = 3 and various potentials. It is apparent that the drivers of 
three-body matrix element evolution depend on the interplay between long- and short-range. 
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FIG. 11: (color online) Contributions from individual terms to the A = 3 ground-state expectation 
value d{vj;^^^ ) / dX for several different two- and three-body potentials, as indicated in the plots. 
We emphasize that A < 2 is very small, comparable to A < 1.5 fm~^ for NN forces in analogous 
calculations with the NCSM 

attractive and repulsive parts. The lower right panel of Fig. [TT] shows an increasing attractive 
strength of the three-body force when starting from an attractive-only two-body potential. 
In this regime the dominant contribution to the evolution of the three-body potential matrix 

(2) (2) 

element is the tree-level two-body connected part [V^ , Vs ]c- This observation accounts 
for the behavior in the right graph of Fig. [9], where the size of the error scales (roughly) like 
(^{2)^2_ Yaj-ying the long-range attraction strength in Fig. [8] shows a similar effect. 

More generally, the impact on ground-state energies of the induced three-body interaction 
depends on the details of the correlations in the wave function. The other plots in Fig. [TT] for 
the more realistic initial two-body potential, V^, show the interplay between two- and three- 
body contributions to the three-body matrix element evolution. The three-body contribution 
to the three-body evolution stays small until the longer-range attractive part of the potential 
begins to affect the evolution. Most of the change is from A = 8 to A = 3, which is 
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dominated by {V , V"*^^^). Thus, the feedback of the three-body potential depends on the 
initial conditions, but in general insures that the binding energy contribution stays small. 

We can repeat the above analysis for A> A and find no fully connected terms with only 
two-body forces. Again, disconnected terms involving two and three body potentials cancel 

(2) (3) 

out in the lower sectors. The leading terms are commutators with one Vr' and one Vs , 
followed by connected terms quadratic in v}^^ and one term quadratic in v}^^ . All terms are 
small and additional cancellations among them further suppress the four-body contribution. 
Thus, the initial hierarchy of many-body forces implies that induced four-body (and higher- 
body) forces will be small. 



IV. SUMMARY 

In this paper we use a one-dimensional system of bosons with short-range repulsion and 
mid-range attraction as a laboratory to explore the evolution of many-body forces with the 
SRG. These calculations serve as a proof-of-principle that induced few-body interactions 
can be calculated without solving T-matrix equations, as necessary in other renormaliza- 
tion group approaches. They establish that working within a harmonic oscillator basis is 
practical, although the generalization to three dimensions will be much more computation- 
ally intensive. This generalization is direct, however, because the recursive construction of 
properly symmetrized few-body basis states used here has already been carried out in the 
context of no-core shell model calculations 0, 0, [HI. The first steps with an A = 3 nuclear 



system will be taken soon, with the resulting three-body matrix elements then used as input 
to calculations in larger nuclei. 

The patterns of SRG evolution observed in our one-dimensional laboratory have many 
similaries to those observed in three dimensions with realistic nuclear interactions. This 
includes, for example, decoupling in few-body systems, which was studied for nuclei with 
two-body forces only [13]. Here we find decoupling is unchanged when induced three-body 
interactions are included. Results for the contributions of induced three-body forces to 
ground-state energies show the same pattern of running observed for nuclear forces when 
a model two-body potential with short-range repulsion and long-range attraction is used. 
Differences in the running when an attractive-only two-body potential is used are understood 
in terms of the dominant contribution to the evolution equation for the three-body potential. 
Starting from an initial decreasing hierarchy of forces (two-body larger than three-body 
and all higher-body forces zero), the hierarchy is maintained for useful ranges of the SRG 
evolution. Induced four- and five-body forces are found to be very small. This is encouraging 
for the applications to realistic nuclei, but of course is not conclusive. 

In a future study, we will use the one-dimensional laboratory to test alternative methods 
to evolve many-body interactions. These include the use of a Slater determinant basis of 
harmonic oscillators ( "m scheme" ) and evolution in momentum space applying the diagram- 
matic methods from Ref. [8|. With the oscillator basis methods we will also consider evolving 
with a choice of the SRG generator Gg that is diagonal in the basis (i.e., Gg = Hose)- The 
block diagonalization characteristic of "V^owfe" potentials (as opposed to band diagonaliza- 
tion as studied here) has been demonstrated for A = 2 in the SRG using a different choice 
for the flow operator Gg- We will test whether this can be extended to A> 3, which would 
circumvent the need to solve coupled T-matrix problems. We will also test an alternative 
to explicit few-body forces, which is to use a normal-ordering prescription in the SRG to 
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generate "density-dependent" two-body interactions. Finally, another avenue to explore is 
the evolution of other operators in few-body systems. 

This work was supported in part by the National Science Foundation under Grant 
No. PHY-0653312 and the UNEDF SciDAC Collaboration under DOE Grant DE-FC02- 
07ER41457. We thank E. Anderson, S. Bogner, R. Perry, L. Platter, J. Drut, and A. 
Schwenk for their comments and discussions. 



APPENDIX A: TRANSFORMATION BRACKETS 

Transformation brackets are the expansion coefficients in the oscillator basis of one system 



of coordinates in terms of another [18|, |19|. We apply them to relate two different choices 
of Jacobi coordinates. Here, we show the relevant transformation using the three-particle 
harmonic oscillator states defined in Eq. and then generalize at the end. 

The single particle momenta are ki, k2, and k^. The unprimed Jacobi momenta [see 
Eq.([3])] are 

Pi = ~ ^2) , 

P2 = ^J^m + k2)/2-k3) , (Al) 
and the primed coordinates are obtained from exchanging ^2 and k^: 

P'l = ^(^^i ~ , 

P2 = l/f ((^1 + h)/2 - k^) . (A2) 

After some algebra, the transformation that exchanges the last two particles (i.e., k2 and 
/cs) can be written as 



Pi 
P'2 




(A3) 



which enables us to express the primed oscillator states in terms of the unprimed ones. 

We denote the three-particle oscillator basis by \n1n2) = Vivll^) where we have set 
hu = 1 for simplicity in this appendix. The transformation that exchanges the last two 
single-particle coordinates can again be written as 



V'l 
V'2 




(A4) 
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The derivation of the harmonic oscillator transformation bracket follows directly as 
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The second line is obtained from operating the transformation on the creation operators 77^. 
The third line is the application of the binomial theorem. The fourth balances the oscillator 
creation and annihilations, and the fifth is just some simplification. 

In the general A-particle system the transformation to exchange the last two particles, 
k^-i and k^, can be written as 



Va^2 
v'a-1 




(A6) 



where d = (A — 1)^ — 1 is the number of generators of the rotation group, U{A — 1), or the 
group U{A) with the center of mass coordinate held fixed. An expression for the bracket 
{nA-2'n'A-iWj^_2n'j^_i) A{A-2), which appears in Eq. is obtained from Eq. flASp by substi- 
tuting the general coordinate transformation Eq. flA6p for the three-particle transformation 



Eq. (jAll), or ^/T/JdTT) and ^/dJJJTd) for 1/2 and ^3/2. 
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